Curvature induced quantum potential on deformed surfaces 
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We investigate the effect of curvature on the behaviour of a quantum particle bound to move on 
a surface. For the Gaussian bump we derive and discuss the quantum potential which results in the 
appearance of a bound state for particles with vanishing angular momentum. The Gaussian bump 
provides a characteristic length for the problem. For completeness we propose an inverse problem 
in differential geometry, i.e. what deformed surfaces produce prescribed curvature induced quantum 
potentials. We solve this inverse problem in the case of rotational surfaces. We also show that there 
exist rotational surfaces in the form of a circular strip around the axis of symmetry which allow 
particles with generic angular momentum to bind. 

PACS numbers: 03.65.-w, 03.65. Ge, 68.65.-k 



The quantum physics of curved surfaces plays an in- 
creasingly significant role in the engineering of devices 
based on coated interfaces pp. Curvature also affects the 
mechanical properties of some biological systems [2.. In- 
vestigation on the geometric interaction between defects 
and curvature in thin layers of superfluids, superconduc- 
tors, and liquid crystals deposited on curved surfaces is 
under way [3j . The curvature of a semiconductor surface 
determines also an interesting mechanism of spin-orbit 
interaction of electrons ]4j. Goldstone and Jaffe[5] and 
Exner[6] proved that at least one bound state exists for all 
two-dimensional channels of constant width except chan- 
nels of constant curvature, which have no bound state. 
The quantum effective potential that appears as a result 
of the presence of curvature in low-dimensional systems 
has been studied also infill EE EH H] and in[T3j where 
it was shown that a charged quantum particle trapped 
in a potential of quantum nature due to bending of an 
elastically deformable thin tube travels without dissipa- 
tion like a soliton. Surprisingly, the twist of a strip plays 
a role of a magnetic field and is responsible for the ap- 
pearance of localized states and an effective transverse 
electric field thus reminisce the quantum Hall effect 14] . 

It is possible to produce very narrow two-dimensional 
conducting surfaces which allow electrons to propagate 
in the channel formed by their boundaries, but require 
the electron wave function to vanish on these boundaries. 
Furthermore, the effects induced by a curved surface on 
the distribution of quantum particles are not fully un- 
derstood. In this paper, we study simple rotationally 
invariant surfaces of varying curvature to gain a broader 
understanding of the interaction between quantum par- 
ticles and curvature and possible physical effects. 
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The results of this paper are based on the exploration 
of the properties of the Schrddinger equation on a sub- 
manifold of M 3 [15j [16] . Following da Costa an effective 
potential appears in the Schrddinger equation which has 
the following form: 

V s (q u q 2 ) = — ^—(M 2 - K )=~(k 1 -k 2 ) 2 , (1) 
2m am 

where m is particle's mass, ft-Plank's constant; q\ and 
q2 are the generalized coordinates on the surface; k\ and 
&2 are the principal curvatures of the surface and M — 
(fej + fc 2 )/2 and K — k\ki are the Mean and the Gauss 
curvatures respectively. 

The presence of the Mean curvature (which cannot be 
obtained from the metric tensor and its derivatives alone) 
in results in an important consequence: V s (qi,q2) 
is not the same for two isometric surfaces, da Costa 
notes [TS]: "thus independent of how small the range of 
values assumed for q% (the third coordinate which mea- 
sures the distance to the surface along the normal vector) , 
the wave function always moves in three-dimensional por- 
tion of space, so that the particle is "aware" of the exter- 
nal properties of the limit surface" . The particle is also 
"aware" of the manner in which it is confined to move 
to that limit surface^ HO]. In view of this, solving the 
inverse problem, we construct a deformed surface which 
corresponds to an effective free one dimensional motion. 
We also report the existence of a circular strip surface 
creating conditions for a zero angular momentum parti- 
cle to bind in a harmonic potential. 

Let us take a rotationally invariant surface r(qi,q2), 
parameterized in Cartesian coordinates in Monge fash- 
ion: 

r(qi,q2) = r(p,(j>) = (p cos (f), p sin <f>, f(p)), (2) 
where p€ [0, oo) and <f> € [0, 2n]. The area element is 

ds 2 = p 2 dcb 2 + (l + /(p) 2 ) dp 2 (3) 
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and the determinant of the metric is given by 

g = p 2 (l + f(p) 2 ), (4) 



where hereafter the dot represents derivative with respect 
to p. From the first and the second fundamental forms 
of this surface the following expressions for the principal 
curvatures are obtained 



Hp) = f(p)(i + f(p) 



-3/2 



Hp) = Hp) p 2 (i + f(p) 



-1/2 



(5) 
(G) 



If the surface is not a plane but it is asymptotically 
planar and cylindrically symmetric then the Schrodingcr 
operator can have at least one isolated eigenvalue of finite 
multiplicity which guarantees the existence of a geomet- 
rically induced bound state|17j. Looking for stationary 
modes in polar coordinates in which the rotational in- 
variance of the surface ([2| is obvious, we separate the 
variables Xt(p, <?M) = exp (— \EtJK) exp (iva.(j))ip(p) (here 
Xt is the tangential to the surface part of the wave 
function [16]) to end up with a quasi-one-dimensional 
Sturm-Liouville equation for the /5-dependent part of the 
wave function: 



<>r i ^=<>r ) IMP) - ^m(p) + ^V'm(p) = ~V s ( P )i> m (p) , (?) 



where V s (p) is given by (flj with ^ and The nor- 
malization condition in p-space is 



2tt 



|^m(p)| yfgdp = 1 



(8) 



with g given by Q. 

Due to the cylindrical symmetry and the conservation 
of the z-componcnt of the angular momentum, we may 
reduce the problem to a one-dimensional equation for 
each angular momentum quantum number m along the 
Euclidean line length along the geodesic on the surface 
with fixed <fi. Introducing the changes [TH] 



^m(p) = F m (x)/y/p, 



(9) 
(10) 



we obtain for the function F m (x) a one-dimensional 
Schrodinger equation which is the Liouville normal form 
of 0: 



dx- 



:F m (x) + [W m (x) - k 2 ] F m (x) = 0. (11) 



Here we have introduced the wave vector instead of the 
energy k 2 = 2mE/h 2 . The geometrical properties of the 
surface determine the quantum effects through the geom- 
etry dependent term W m (x) in equation (UTJ): 



(12) 



W m [x(p)]^~ 1 -k 2 1 (p)+ m2 2 1/4 , 
4 p z 



where k\ is given by The normalization condition in 
x-space is 



/>oo 

2tt / \F m (x)\ 2 dx = 1. 
Jo 



(13) 



The term in W m , proportional to m 2 , is the potential 
that describes the familiar centrifugal force. Less famil- 
iar is the negative correction term —1/4 which results 
not from the angular motion but from the radial motion 
(and can be traced back to the radial derivatives in the 
Laplacian expressed in the associated with the surface 
coordinates ([2])). This is a coordinate force that comes 
from the reduction of space from three to two dimensions 
and was called quantum anti- centrifugal force by Cirone 
et. al. |21j because it possesses binding power due to 
quantum mechanics. Similar situation for the free radial 
motion of a particle was noticed in|22j. This contribu- 
tion is strengthened by the binding curvature induced 
potential — fc 2 /4, a geometric force. 

An interesting situation (that may have practical ap- 
plication) unfolds when we consider a Gaussian bump 
which appears often when considering deformations of a 
surface with the following profile: 



f(p) = -A CXp (-p 2 /cr 2 ) 



(14) 



where do is the dispersion of the profile and Aq > is its 
depth. Thus prior to the attempt to search for a solution 
of (111 with (12) we would like to discuss the physics. 

The potential W m is repulsive for particles with non- 
vanishing angular momentum (m ^ 0) and no solutions 
with negative energy exist. 

The effect from different contributions in the potential 
W m stands out most clearly for particles with zero angu- 
lar momentum, that is, m = 0. These are attracted to 
the origin and are found in a ring-shaped region around 
the axis of symmetry, while all particles with m^O are 
repelled from the center. It is clear that the bump not 
only introduces scale but also breaks the symmetry of 
the R 2 plane introducing a natural origin - its center. 

The Gaussian bump (14 1 also introduces a character- 
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FIG. 1: A stretched strip on a bump and its downward pro- 
jection on a flat plane. 



istic length on the surface. By vertically stretching the 
surface, the amount of area in the stretched strip is in- 
creased with respect to its projection on the plane which 
is depicted in Figure ([TJ. As a result, the quantum par- 
ticle would rather remain localized when it is situated 
near the maximal stretch of the surface thus introduc- 
ing a scale of the order of the dispersion a of the pro- 
file (14). This can easily be seen when looking at the 

h: 5x is bigger 



Heisenberg uncertainty relation 5p5x 
in the stretched area. Hence, it is energetically favorable 
the quantum particle migrates to the stretched area of 
the bump thus creating charge distribution if it itself is 
charged. This charge density is exactly depicted by the 
probability density of finding the particle in a specific 
place on the surface. 

Since we are interested in bound states now we con- 
sider only particles with m = 0, since only for them 
the potential Wq is negative and they will bind with dif- 
ferent quantized energies. Equating the kinetic and po- 
tential energy terms we provide a rough estimate of the 
order of magnitude of the bound state energies, namely 
E < — ti 2 /2mo-Q. Later on in the text we will see that 
this estimate is correct up to a geometry dependent fac- 
tor Aq/ctq, which measures the flatness of the profile. 



Obtaining a solution to equation (11) is not an eas- 
ier task than solving the original one. Nevertheless, we 
can be sure that bound states do exist since we have 
transformed the problem to a one-dimensional one with 
a negative potential [1 9 . Using an approximation we can 
transform this equation into simpler one thus providing 
us with a chance to generate an approximate wave func- 
tion. The behaviour of the wave function in the vicinity 
of the origin, that is p — » can be considered as a start- 
ing point in our consideration. In that particular limit 



p due to rt9| and lim p _^ / = (for the Gaussian 



bump) for the potential Wq(x) we have 



Wo[x] 



1 

Ax 2 



3± 

7 2 ' 

CT o CT o 



(15) 



Equation (111 with the above potential possesses a stable 
negative energy solution with k' 2 = 2m\E\/h 2 — A^/a^ 
provided 2m\E\/h 2 > Al/a^. Here A^/a 2 is a measure 
of the flatness of the Gaussian bump. Due to the sign 
change of the energy the ordinary Bessel functions Jo 
and Yo which solve equation (11) with the above poten- 
tial ( |15| turn into the modified Bessel functions Jo and 
Kq 20 . Since the modified Bessel function I increases 
exponentially for large distances, whereas Kq decreases, 
the boundary conditions imposed by the asymptotic flat- 
ness of the Gaussian bump and the need for a squared 
integrable wave function, we select K . Thus the solution 
in the vicinity of x — + is expressed with the modified 
Bessel function 

F {x) sa y/x~K (k'x). 

Having deduced the correct behaviour of the wave 
function in the vicinity of the origin we can try to solve 
equation (11) with the ansatz 



Fq(x) = cxp 



:S(X) 



fxK {k'x). 



(16) 



Using the fact that Kq satisfies the Bessel equation we 
end up with an approximate equation for S{x) 



-S 



h (4 



h x 



J2 
A 



(17) 



where hereafter ( ) denotes differentiation with respect 
to a;. In the derivation of ( 17 ) we kept in mind that x sa p 



and approximated 1/p 2 with 1/x 2 in the potential ( |12[ ). 
We have also approximated dKo/dx w —Kq/x which 
preserved the behaviour near the origin. 

Next we expand S(x) in series with increasing powers 
of the small parameter h a la WKB type expansion 

S(x) = Sq(x) + ihS\{x) + . . . 

and obtain for the first two functions in the expansion 
the equations 



S = => So = const, 



S, = (Srf + SJx- 



J2 

A) 



V 



(18) 



(19) 



Solving (19 1 we assume that i.) 
(Si) 2 > S 1 /x to obtain 



(S x ) 2 > Si and ii.) 



S 1 = ± 



Al 



1 

4 



\°* 4 (l + / 2 )' 



(20) 
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From the form of the solution we see that the assumptions 
i.) and ii.) we made in deriving it are justified away from 
the origin. Since we already dealt with the behaviour of 
the wave function in the vicinity of that point we write 



F (x) 



e iS /h 

V2ttN 2 



\fxKo(k' x) exp 



Si dx 



(21) 



where we have chosen the + sign solution of ( 20 I , ame 



liorating the vanishing behaviour of the wa ve f unction at 
infinity. The norm N 2 is calculated using (13). The log- 



arithmic divergence at the origin of the modified Bessel 
function of second kind Kq is regularized by the extra 
factor proportional to y/x in terms of the Euclidean line 
length on the surface, making the probability density to 
find the particle at the origin vanishing. Moreover, since 
(21 1 vanishes exponentially for large distances, the ra- 
dial probability displays a maximum close to the origin, 
which is exactly the behaviour of the expectation value 
that we already deduced from the considerations using 
Hcisenberg's principle. 

As a result, particles with m = concentrate in a 
ring shaped region on the surface of the Gaussian bump 
around the axis of symmetry. 




FIG. 2: The probability density \4>o(p)\ 2 y/g (vertical axis) to 
find a particle between p and p + dp (horizontal axis) for the 
Gaussian bump. Here we have set Ao = <jo = 1 and k' = 5. 
The initial phase So = 0. 

Figure Q depicts thus generated probability density 
\^o{p)\ 2 y/g of finding the particle in a state with angular 
momentum m = between p and p + dp in terms of the 
distance from the axis of symmetry and clearly shows 
that the probability of finding the quantum particle is 
biggest at the stretched area of the bump as already dis- 
cussed. This shows that the approximations which lead 
to (21) are fully justified. They also show the way to a 



proper test function for a run of the variational method 
to determine the ground state energy and wave function, 
which must not have nodes - a property clearly visible 



in (21 ). This is a matter of a future study. 



Here we do not present an estimate of the structure 
of the energy spectrum except the already stated fact - 
the discrete spectrum starts at E < —f?A\l2ma\, value 



which is determined by the geometrical extensions of the 
Gaussian bump. Over that value the motion of a particle 
with m = is classically infinite and hence the wave 
function is supposed to be in the form of a plane wave. 
What is, at least, the approximate expression of that 
plane wave? 



We propose to look for solutions of (11) in the form of 



K T=/o(z)exp ±t / Q dx 



(22) 



where Q w y/kj/A + 1/Ax 2 + k 2 . Here f (x) 
is fast changing function, such that F (it) 



y—Qo ± 2jQ /o//o ± iQo ) ^b- From here and (|ll|) we 

see that 



fo ~ 

Hence we write 



lQo 
2Q 



fo ~ lh 



TTicont 
^0 



(X) 



exp 



(±i J y/kf/A+l/Ax 2 + k 2 dx^ 
(k 2 /A+l/Ax 2 + fc 2 ) 1/4 



(23) 



(24) 



where the plane wave exhibits a node in the origin. 

The energy spectrum for particles with non-vanishing 
angular momentum is continuous and starts at E = 0. 
Good approximate solutions of (111 are the ordinary 



Bessel functions y/xJ m (kx) and y/xY m (kx), whose ar- 
gument is the Euclidean line length on the surface. 

Let us now turn our attention to the inverse problem 
or equivalently the question "What rotationally invariant 
surface leads to an effective geometry induced potential 
W m that equals prescribed negative function —U[x(p)], 
where U > for Vp?" (negative because we are primary 
interested in bound states). This question is of particu- 
lar interest since for certain classes of negative potentials 
— U we already know the exact wave functions which can 
readily be used in revealing the particle's distribution on 
the surface. The solution of the inverse problem goes 
through the recognition of its equivalence with the fol- 
lowing differential equation (see equations (12 1 and (pi)): 



hp? 



i + f(py 



U{ P ) 



1/4 



(25) 



From equations ^ and (25) one can easily deduce a 
condition on k\(p) in order to have a binding potential 
for a particle with m/0, i.e. 



k\ > 



4nr 



For smooth surfaces at p = this means that there will be 
only strips where this condition holds. For flat surfaces 
where ki = this condition implies that only particles 
with m = bind, a fact previously noticed inj^U 122"] . 
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Now we try to solve equation (25 1 using substitution 



/ = sinh (oj) (here to — uj(p)) to end up with the lin- 
ear equation 



1 dtanh(w) 

2 dp 



±JU(p) 



1/4 



(26) 



which can be solved thus yielding a result for f(p). 
integrate to obtain the profile of the surface 



We 



f(p) = ± f 
J pi 



\A\ 
\/l~42 



dp', 



(27) 



where 



a( P ) = ±2 r 

■J Po 



lU(p') 



m- 



-1/4 



dp' . (28) 



Here po and pi are constants of integration and are to 
be determined by the boundary conditions due to the 
behavior of the function U representing the potential we 
want to model. Since f(p) takes only real values (the 
same is true for all of its derivatives, i.e. d n f/dp n E R 
for n = 0, 1, . . .) as a function describing the profile of 
a surface we impose < |„4| < 1. The A = Q case is 
realized by a flat surface. Using the theorem of the mean 
value in ( 28 1 and the above inequality we obtain 



po < p < po + 



1 



u(0 



m 2 - 1/4 



-1/2 



(29) 



where the point £ £ [po , p] . In that manner we show that 
for generic angular momentum m and non-zero poten- 
tial U the corresponding rotational surface creating this 
potential and allowing a bound state of a quantum par- 
ticle with m ^ exists only in a ribbon, a circular strip 
around the axis of symmetry. 

Let us also note that for m ^ the bump (a macro- 
scopic structure) has a magnetic moment, i.e. the macro- 
scopic deformation of the surface acquires quantum num- 
ber. Indeed the probability density current J (div J = 0) 
associated with the wave function xt is given by 




lie fegg3k ,0 

1 + /2 



where J^, J p and J z are the components of the vector 
J in the orthonormal right-handed triad (e p ,e^,e z ), has 
nonzero and quantized with m circulation along the cir- 
cumference of the bump. 

Now let us give a couple of examples: 



• U — 0, i.e. free motion; From (28 1 and (27 1 it 
follows that 



Jpi/po J i - (4m 2 - 1) In 2 (p 1 ) 



(30) 




FIG. 3: The surface f(p)/po (vertical axis) with a cusp on 
which a particle with m = moves freely. On the horizontal 
axis is plotted p/ po- At infinity that surface tends to a cone. 



Here po determines the characteristic scale of the 
surface on which we consider this free quantum 
problem. Its magnitude can be chosen in [pm] or 
[nm] depending on the scale that we want to model. 
If we consider the shape of a surface allowing a free 
motion of a particle with zero angular momentum, 
that is m = 0, we can integrate ( |30| with p\ = po to 
produce Figure pj. That surface is asymptotically 
flat as it tends to a cone fo cc (p) — > p as p — > oo. 

Next we consider the strip of a surface allowing a 
free motion of a particle with non-zero angular mo- 
mentum, that is m 7^ 0. We can impose the condi- 
tion 1 > (4m 2 — 1) In 2 {p/ po) to find the extensions 
of the strip 



po < p < po exp (4m 2 - 1) 



-1/2 



(31) 



Since the problem is defined on a strip, we can 
quantize it in the usual finite volume method, which 
would lead to standing wave solutions on the sur- 
face. Their energy is 



mpl [t 



(4m 2 



in 

-l)-!/2 



-1 



(32) 



• U = uJ 2 p 2 , i.e. harmonic oscillator potential; From 
(|28| it follows that 



A h ^(p) = ±2l ^V- + 

'pa 



P" 



(33) 



Let us consider two cases consequently: 

I.) The m = case yields for (33 1 the result 



arctan 



(4w 2 p 4 - 1) 



-1/2' 
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where we have set the value of po = l/y/2\u>\. From 
the requirement < |.4.Q arm | < 1 we obtain the 
following estimate 



(l + (2 + e) 2 ) 1/4 _ 5V4 + E/5 3 / 4 



P < 



(35) 



where e = ir/2 — arctan (4w 2 p 4 — l) "<1. Thus 
we find an expression for the extensions of the cir- 



cular strip of a rotational surface creating harmonic 
potential and allowing a bound state of a particle 
with vanishing angular momentum in a harmonic 
oscillator potential 



1 



<P< 




(36) 



II.) The m ^ case yields for (33) the result 



± A^ m (p) = V^V + m 2 - 1/4 - Vm 2 - 1/4 arctanh 



m 2 - 1 /4 



-1/2' 



(37) 



r 



Next we impose the condition < |-4^ rm | < 1 and 
numerically find for m = 1, 2, 3 particles 



1.071 < V2uj P < 1.602 
1.602 < y/2up < 1.957 
1.980 < V2ujp < 2.265 



m = 1, 
m = 2, 
m = 3. 



(38) 



In conclusion we have calculated the quantum poten- 
tial associated with a Gaussian deformation of a two- 
dimensional plane. For that simple surface deformation 
we have shown that a quantum binding force of geometric 



origin attracts only particles with m = to the origin. 
Further, we have solved the inverse problem and have 
shown that surfaces in a form of a ribbon allow bound 
states of particles with generic angular momenta. Thus 
we speculate that a classical object (the ribbon) exhibits 
quantum characteristics (the magnetic moment due to 
non-vanishing quantized probability current circulation 
along the circumference) acquired due to curvature. We 
have also found and depicted the surface which corre- 
sponds to a free motion of a particle with vanishing an- 
gular momentum. 
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